Universal behavior of pair correlations in a strongly interacting Fermi gas 



E. D. Kuhnle, H. Hu, X.-J. Liu, P. Dyke, M. Mark, P. D. Drummond, P. Hannaford, and C. J. Vale 

ARC Centre of Excellence for Quantum-Atom Optics, 
Centre for Atom Optics and Ultrafast Spectroscopy, 
Swinburne University of Technology, Melbourne 3122, Australia 
(Dated: June 29, 2010) 

We show that short-range pair correlations in a strongly interacting Fermi gas follow a simple 
universal law described by Tan's relations. This is achieved through measurements of the static 
structure factor which displays a universal scaling proportional to the ratio of Tan's contact to the 
momentum C/q. Bragg spectroscopy of ultracold 6 Li atoms from a periodic optical potential is used 
to measure the structure factor for a wide range of momenta and interaction strengths, providing 
broad confirmation of this universal law. We calibrate our Bragg spectra using the /-sum rule, 
which is found to improve the accuracy of the structure factor measurement. 



PACS numbers: 03.75.Hh, 03.75.Ss, 05.30.Fk 

Universality is a remarkable property of strongly in- 
teracting systems of fermions [IH5] ■ Universality means 
that all dilute Fermi systems with sufficiently strong in- 
teractions behave identically on a scale given by the aver- 
age particle separation. With the discovery of universal- 
ity in the Bose-Einstein condensate (BEC) to Bardeen- 
Cooper-Schrieffer (BCS) crossover, ultracold Fermi gases 
near Feshbach resonances have become a central topic in 
atomic physics |¥HTU]. Strongly interacting Fermi sys- 
tems arise in a wide variety of settings, from astrophys- 
ical to nuclear and condensed matter systems. One can 
therefore study universality in ultracold atomic gases to 
help understand other strongly interacting Fermi super- 
fluids, taking advantage of the ability to precisely control 
the atom-atom interactions. 

Understanding these strongly interacting Fermi gases, 
however, poses significant challenges [11 . In 2005 Shina 
Tan made dramatic progress by deriving several exact re- 
lations for Fermi gases in the BEC-BCS crossover, which 
relate the microscopic properties to bulk thermodynamic 
quantities |12ffl4] . These exact relations are applicable in 
broad circumstances: zero or finite temperatures, super- 
fluid or normal phases, homogeneous or trapped systems, 
and in few or many-body systems. 

In this letter, we experimentally verify a universal re- 
lation for short-range pair correlations |12j using Bragg 
spectroscopy. This is achieved through measurements of 
the static structure factor, given by the Fourier trans- 
form of the pair correlation function |15| . The structure 
factor of a unitary Fermi gas has an exact scaling with 
the ratio of Tan's contact to the momentum C/q. For 
systems with finite scattering length, we also confirm the 
first order correction to the universal law. Our measure- 
ments are compared to new calculations for the contact 
based on a recently developed below-threshold Gaussian 
pair fluctuation theory |16] . 

The contact C in a two component Fermi gas quanti- 
fies the likelihood of finding two fermions with opposite 
spin close enough to interact with each other. In systems 



where the range of the interaction potential is negligible, 
this single parameter encapsulates all of the information 
required to determine the many-body properties [TTl ITS] . 
C depends on the s-wave scattering length, density and 
temperature of the system. Tan showed that the inter- 
nal energy of a gas across the BEC-BCS crossover can 
be expressed functional of the momentum distri- 

bution which has a C/g 4 dependence at large momen- 
tum q and that the pair correlation function diverges as 
C/r 2 at short distance r < 1/fcf, where Uf is the Fermi 
wavevector |121 ITS] . Tan also derived the adiabatic rela- 
tion dE/d(—\/a) = h 2 C/(4:TTm), where m is the atomic 
mass, giving the change in the total energy E due to 
an adiabatic change in the scattering length |13| and ex- 
tended the virial theorem to finite a and imbalanced mix- 
tures [2] . The contact C was first extracted [TH| from the 
number of closed-channel molecules determined through 
photo-association [5] and the adiabatic and virial Tan 
relations were very recently verified experimentally |20| . 
We will generally refer to the dimensionless contact X 
given by C/(Nkp) where N is the number of particles. 

Short-range structure in a quantum fluid depends upon 
the relative wave-function of the interacting particles, 
in this case fermions in different spin states. In a two- 
component (spin-up/spin-down) Fermi gas with contact 
interactions this is given by ?/>^(r) cx 1/r — 1/a, where a 
is the s-wave scattering length. Starting from this wave- 
function, Tan showed that the spin-antiparallel pair cor- 
relation function is given by Eq. (1) which includes the 
contact as a pre- factor [T^] 

Pair correlation functions are difficult to measure di- 
rectly in ultracold gases; however, it is possible to mea- 
sure macroscopic quantities which depend on correlation 
functions in a well defined way. A prime example is the 
static structure factor, S(k), which is given by the Fourier 
transform of (r) (q = hk is the probe momentum) . In 



2 



a two component Fermi gas with an equal number of par- 
ticles in each state N, the structure factor consists of two 
components, corresponding to correlations between par- 
ticles in the same state and particles in different states 
S(k) = S<ft(k) + Sfi(k). When the momentum is much 
larger than the Fermi momentum, particles in the same 
state will be uncorrelated and Sff(k 3> kp) ~ 1 so all 
variation in S(k) will then be due to changes in S^i(k) 
[15, 29 . The Fourier transform of Eq. 1 yields the follow- 
ing for S n (k) 
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which has a straightforward dependence upon Z, the rela- 
tive probe momentum k/kp, and the interaction strength 
l/(kpa). At unitarity, a — > oo, the second term vanishes 
and S(k) varies linearly with kp/k. Rewriting Tan's re- 
lation in this way points to a method to experimentally 
verify Tan's universal relation for pair correlations be- 
tween spin-up/spin-down fermions since the static struc- 
ture factor can readily be measured using Bragg spec- 
troscopy. 

Apart from the contact, all parameters in Eq. (2) are 
easily determined or set by experimental parameters. We 
employ new methods to calculate the contact which over- 
come the need for interpolation schemes between the lim- 
iting BEC and BCS cases used in earlier studies [T3]. To 
find the low temperature contact we first calculate the 
ground state energy across the BEC-BCS crossover using 
a recently developed below-threshold Gaussian pair fluc- 
tuation theory [16_, which has shown excellent agreement 
with thermodynamic measurements [9, 10, 21 and quan- 
tum Monte Carlo (QMC) simulations [25]. We then use 
Tan's adiabatic theorem to find the contact. The high 
temperature (T > 0.5 Tp) contact is calculated from a 
quantum cluster expansion |23| . Figure [I] shows the cal- 
culated contact X for a range of temperatures through the 
BEC-BCS crossover, along with the results from photo- 
association (Fig.[lji) [HMTH]- The contact increases mono- 
tonically from the BCS to BEC regimes and is highest 
at low temperatures. Combining these calculated results 
for the contact with Eq. (2) for high momentum transfer 
Bragg scattering, we obtain a direct, quantitative predic- 
tion that is readily tested experimentally. 

Inelastic scattering experiments are a well established 
technique to probe both the dynamic and static structure 
factors of many-body quantum systems |24| . Ultracold 
atoms are highly amenable to inelastic scattering through 
Bragg spectrosocpy which has previously been used to 
measure both the dynamic |2"5l |2"6"| and static structure 
factors |27l |2"5] of atomic Bose-Einstein condensates. In 
ultracold Fermi gases Bragg spectroscopy was used to 
measure both the dynamic and static structure factors 
over the BEC-BCS crossover, albeit at a single momen- 
tum |29j . The ability to vary the momentum differen- 
tiates Bragg spectroscopy from rf spectroscopy [3TJI 1ST] , 
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Figure 1. Theoretical prediction for Tan's contact in the BEC- 
BCS crossover, (a) Evolution of contact for a trapped Fermi 
gas with l/(fci?a) at different temperatures. Points are exper- 
imental data from refs [51 119|. The T = curve is calculated 
using the contact for a homogeneous Fermi gas (shown in b) 
and a local density approximation. Finite temperature con- 
tact is determined from a cluster expansion theory |23| . (b) 
Zero temperature contact in free space, obtained from a Gaus- 
sian pair fluctuation theory |16| . Red circle shows the contact 
from QMC for pair correlation functions at unitarity |15j . (c) 
Temperature dependence of the trapped contact at unitarity. 
With decreasing T, contact increases rapidly towards the T 
= result (purple star). 



providing access to a broad range of the excitation spec- 
trum and avoiding final state interaction effects as no 
third atomic state is involved. 

To demonstrate the validity of Eq. (2) we need to mea- 
sure S(k) for a range of momenta (k/kp) and at different 
values of the interaction parameter l/(kpa). The start- 
ing point of our experiments are clouds containing N ~ 
3 x 10 5 6 Li atoms in an equal mixture of the two lowest 
lying ground states \F — 1/2, rap — ±1/2) evaporatively 
cooled in a single beam optical dipole trap (A = 1075 
nm) at a magnetic field of 834 G [32]. Next, we adia- 
batically ramp up a second far detuned laser (A = 1064 
nm) which intersects the first trapping beam at an angle 
of 74° forming a crossed beam dipole trap. By appro- 
priately selecting the intensities of each of the two trap 
lasers we tune the mean harmonic confinement frequency 
of the crossed trap u> over the range uj = 2ir x (38 — > 252) 
s _1 while the aspect ratio, Lj xy /oj z , varies from 3.4 to 16. 
Controlling the trap frequencies in this way allows us to 
tune the atom density and hence the Fermi wavevector 
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over a broad range. Once in the final crossed beam trap, 
we ramp the magnetic field to select the scattering length 
a which gives the desired value of l/(kpa). This allows us 
to tune kp and l/(kpa) independently. The cloud is then 
held at the final magnetic field for a time r 10w _1 to 
reach equilibrium before applying the Bragg pulse and 
imaging at that magnetic field. 

Bragg scattering is achieved by illuminating these 6 Li 
clouds with two counter propagating laser beams with 
a very small frequency difference S. This creates a 
standing wave which moves with velocity 5/ks r where 
ksr — 27r /Xsr is the wave- vector and As r = 671 nm 
is the wavelength of the Bragg lasers. When the en- 
ergy difference between the two Bragg lasers equals the 
kinetic energy associated with a two-photon recoil, res- 
onant Bragg scattering occurs (i.e. S = 2hk 2 Br / m) . The 
Bragg lasers are approximately 2 GHz red-detuned from 
the two ground hyperfine states which is large compared 
to the ~ 80 MHz splitting between the 1 1/2, ±1/2) states. 
This means that both states are coupled approximately 
equally to the Bragg beams (to within 4%). 

To measure the static structure factor S(k) we begin by 
recording a Bragg spectrum. This consists of measuring 
the momentum transferred to the cloud by the Bragg 
lasers as a function of 5. The transferred momentum 
is directly proportional to the resultant center of mass 
displacement AX(k, 5) measured after 2 ms time of flight 
[2"9"] . A Bragg spectrum for a unitary gas is shown in the 
inset of Fig. 2 for k/kp — 8.5 as a function of the Bragg 
frequency S. At this momentum, the pair and free-atom 
excitations are clearly distinguished at frequencies of ~ 
150 kHz and ~ 300 kHz, respectively. The Bragg pulse 
duration, 50 ^s, is short compared to the two-photon 
Rabi cycling period, ensuring that spectra are obtained 
in the linear response regime. 

The center of mass displacement, AX(k, S), is propor- 
tional to the convolution of the spectral content of the 
Bragg pulse and the dynamic structure factor of the gas 
S(k,5) [33 . The proportionality constant depends upon 
the two-photon Rabi frequency flsr which can be dif- 
ficult to accurately measure. Integrating AX(q, 8) over 
d yields a number proportional to the static structure 
factor, S(k) = HN~ l J S(k,S)d8, but £lsr remains un- 
known. We overcome the need to find Qsr by invok- 
ing the /-sum rule (22]: NE r — J S(k,6)8dS, where 
E r = 2h 2 k 2 Br /m is the two-photon recoil energy. As both 
S(k) and the /-sum rule involve N, we normalise the area 
under each measured spectra using 



S(k) = 



2hk 2 Br fAX(k,5)dS 
~^jAXlj^)6d5 : 



(3) 



where the constant involving Qsr appears before each 
integral and therefore cancels. Equation (3) is absolute, 
requiring only knowledge of the recoil energy, which can 
be determined with high accuracy, leading to an accurate 
measure of S(k). In Fig. 2 we plot S(k) obtained in this 
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Figure 2. Static structure factor in the BEC-BCS crossover. 
Experimental points were obtained by integrating Bragg spec- 
tra normalised via the /-sum rule for k/kF = 4.8. Error bars 
are due to shot to shot atom number fluctuations and un- 
certainties in measuring center of mass. The theoretical line 
is a zero temperature result, calculated by interpolating the 
near-resonance structure factor Eq. (2) with the asymptotic 
results in the BCS and BEC limits. Inset: A Bragg spectrum 
obtained at l/(kpa) = and k/kF = 8.5 showing center of 
mass displacement AX(k, 8) versus Bragg frequency S. Points 
are single shots and the line is a guide to the eye. The pair 
(5 ~ 150 kHz) and free atom peaks (5 ~ 300 kHz) are clearly 
distinguished. 



way for k/kp = 4.8 which shows excellent agreement with 
the zero temperature theoretical calculation. 

Next we come to the measurement of S(k) as a func- 
tion of k/kp. At unitarity we expect S(k) to vary linearly 
with kp/k according to Eq. (2). Rather than changing 
k directly by varying the angle between the two Bragg 
beams, we change kp = (48A r ) 1 / 6 \/mw /h by varying 
u) = (cj^WyWj) 1 / 3 , the geometric mean frequency of the 
optical dipole trap. Universality allows us to change the 
relative length scale being probed simply by changing 
the density of the gas. Using the crossed beam opti- 
cal trap described above we can tune w such that the 
Fermi wavevector can be tuned anywhere over the range 
&f = 2.1 /xm _1 — > 5.3 iMn^ 1 , or k/kp = 3.5 — > 9.1. 
A sequence of Bragg spectra were obtained for gases 
at three values of \/(kpa) = +0.3, 0.0 and -0.2, while 
varying kp. From these, S(k) was extracted and the 
results are plotted in Fig. 3. The solid lines are the pre- 
diction from Eq. (2) using the zero temperature contact 
with no free parameters. The experimental points closely 
follow the theory, confirming the exact analytic depen- 
dence of S^i(k) = S(k) — 1 on q, Eq. (2). At unitarity 
1 / (kpa) = 0, the dependence on a vanishes and a straight 
line fit (dashed line) shows the simple universal behav- 
ior of fermionic pairing. The fitted slope of 0.75 ± 0.03 
is slightly below the zero temperature prediction of 0.81 
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Figure 3. Universal dependence of the static structure fac- 
tor of a strongly interacting Fermi superfluid. Measured and 
calculated static structure factor versus fcf/fc for l/(fcj?a) = 
+0.3, 0.0, and -0.2. Bragg momentum k is fixed while &f is 
varied by changing the mean trapping frequency Q. Vertical 
error bars are due to atom number fluctuations and uncer- 
tainties in measuring the center of mass and horizontal error 
bars are due to atom number fluctuations and uncertainties in 
Q. Solid lines are the zero temperature theory and the dashed 
line is a straight line fit to the l/(kpa) = data yielding a 
slope of 0.75 ± 0.03. 

which may be due to reduced pairing at the finite tem- 
perature (T/T F = 0.10 ±0.02 at unitarity). At l/{k F a) 
= -0.2 the temperature will be lower following the adia- 
batic magnetic field sweep [34] while at \/(kpa) = +0.3 
the temperature will be higher but pairing takes place 
at much higher temperatures. At T <C Tc (— 0.2 Tp), 
phonons dominate the excitations and the contact should 
increase as (T/Tp) 4 [35]. In the relevant temperature 
window, we estimate this increase to be only 0.1% which 
could easily be negated by normal single-particle excita- 
tions localized at the cloud edge. At l/(k F a) = +0.3 the 
data depart from a straight line displaying the downward 
curvature consistent with the first order term in Eq. (2). 
A similar upward curvature is seen at l/{kpa) = -0.2. 
Our simple relation Eq. (2) is seen to accurately describe 
S(k) on both sides of the Feshbach resonance demon- 
strating the wide applicability of the Tan relations. 

In summary, we have shown that the structure fac- 
tor of a strongly interacting ultracold Fermi gas follows 
a universal law which is a direct consequence of Tan's 
relation for the pair correlation functions. Our measure- 
ments provide one of the first demonstrations of a broadly 
applicable exact result for Fermi gases in the BEC-BCS 
crossover. This work opens the way to a complete tem- 
perature and interaction dependent map of the contact 
through the BEC-BCS crossover and may provide a new 
means for obtaining the equation of state. 
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